Heavy quarkonia description from an energy dependent quark-antiquark potential. 
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We model the spectral effect of open flavor meson-meson thresholds in heavy quarkonia. The 
proposed energy dependent quark-antiquark static potential tries to incorporate in a quark model 
scheme the results from unquenched lattice calculations. A good qualitative and a reasonable 
quantitative description of charmonium and bottomonium, including the "new charmonium-like" 
states, is obtained. 
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I. INTRODUCTION 

The discovery of A (3872) [3-[| opened a new era 
in heavy quarkonia spectroscopy since its properties 
could not be reasonably reproduced with any conven- 
tional cc quark model employed until then for the de- 
scription of charmonium. Since its mass is about the 
D°(1865)D (2007)° threshold energy alternative descrip- 
tions based on tetraquark or meson-meson models have 
been proposed to explain its nature. In this regard it 
should be pointed out that the J PC quantum number as- 
signment to A (3872) is still an experimental issue. Other 
resonances such as A (4260), A" (4360) ... that were dis- 
covered later are considered unconventional cc states as 
well (see and references therein). 

Taking for granted that these resonances are not pure 
conventional cc states and realizing that all of them are 
close to open flavor meson-meson channels it is impor- 
tant to analyze the effect of these meson-meson thresh- 
olds on the quark-antiquark interaction and the role they 
could play in the description of such non conventional 
states. This effect has been studied in a general form 
in lattice QCD by considering the static quark-antiquark 
(QQ) ground state energy when an open flavor meson- 
meson system may be induced by sea pairs (qq). The 
results are summarized in Fig. 22 of reference [7j that we 
reproduce here for completeness as Fig. [TJ 

The two thin curved lines following lattice data (circles 
and pentagons respectively) represent the adiabatic sur- 
face mixing produced by the presence of one static-light 
two meson threshold gcncrically called BB whereas the 
three thick lines correspond to an educated guess for the 
case of two thresholds generically called BB and B S B S . 
When not close to threshold the quark-antiquark ground 
state energy E{r) (r is the quark-antiquark distance) cor- 
responds to the three-parameter fit or — - + -Eq with the 
values of the parameters differing in the one and two 
threshold cases. In particular the string tension a di- 
minishes when going from one to two thresholds. When 
close to a threshold the expression for E{r) gets modi- 
fied due to string breaking: the corresponding eigenstates 
contain threshold state components (BB, B S B S ...), the 
closer the energy to that of the threshold the higher its 
probability. 

We shall use these results from lattice to build in Sec- 




FIG. 1: Calculated energy levels from lattice QCD for two 
sea quark flavors: circles and pentagons over the thin lines. 
Speculation for three sea quark flavors: thick lines. 



tion II a non relativistic quark model based on an energy 
dependent quark-antiquark static potential V(r). As this 
model will extend the conventional quark model to incor- 
porate effects from meson-meson (threshold) components 
we shall call it Extended Quark Model (EQM). Then in 
Sections III and IV we shall apply it to the description 
of charmonium and bottomonium respectively. Finally 
in Section V our main results and conclusions will be 
summarized. 



II. EXTENDED QUARK MODEL (EQM) 

Aiming at a description of the heavy quarkonia (Q = 
b, c) spectra and properties we shall rely on a non rela- 
tivistic quark potential model. We expect this to be valid 
for bottomonium and to a certain extent (up to relativis- 
tic corrections) also for charmonium. Heavy quarkonia 
masses will be identified with the energy eigenvalues of 
the QQ system obtained by solving the Schodinger equa- 
tion with an energy dependent quark-antiquark poten- 
tial. We shall assume that the effectiveness of the pa- 
rameters (quark masses and parameters of the potential 
fitted to reproduce the spectra) may be appropriately 
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taking into account, at least in part, kinetic and poten- 
tial energy corrections. 



A. Energy Dependent Potential 

In the non relativistic static approximation the quark- 
antiquark ground state energy E(r) may be identified 
with the sum of the masses of the quark (tuq) and the 
antiquark (tu-q) plus the static quark- antiquark poten- 
tial V(r). Let us consider the more general two thresh- 
old case (thick lines in Fig. [1]). The resulting potential 
is energy dependent because its form, dictated by E(r), 
varies according to the spectral region considered. 

Thus for quark-antiquark spectral energies below the 
mass of the first threshold E spec < (M^ — 2ms 

in Fig. [1]), the potential V(r) has at small distances the 
form (or — ~ + Vo). We shall include the constant V in 
the definition of the quark and antiquark masses so that 
we shall write the potential as or — -. This form main- 
tains up to a distance close below the crossing distance 
tt % defined from Vlr^ ) = ott x — = Mt x — inQ — rn-g. 
Then V(r) starts to flatten approaching its asymptotic 
value in this energy region M<r x — tuq — m-Q. The phys- 
ical mechanism underlying this flattening of the quark- 
antiquark interaction is the creation of a qiqi pair that 
recombines with the static quarks (QQ) giving rise to the 
two static-light mesons OTiTj ((Qqij) an d 2^2Ti ((Qqi)) 
so that Mt ± = M<rjyi 1Tl + M<xn 2Ti . We shall call this two 
meson threshold state |Ti) . 

For the spectral region Mt x < E spec < Mt 2 where 
Mt 2 stands for the mass of the second threshold (2m b s 
in Fig. [1]) V (r) is equal to — ttlq — m-Q- up to a dis- 
tance close below , then it rises until getting for a 
distance close above tt x the form or — -. This form 
is maintained up to a distance close below the cross- 
ing distance ry, defined from V(rT 2 ) = or^ 2 — = 

~ " r T 2 

Mt 2 — rag — m-Q where V(r) starts to flatten approach- 
ing its asymptotic value Mt 2 — — tUq. The underly- 
ing physical mechanism implies the annihilation of q{q~i 
from (Qq~i) miTl ~ (Qqi) m ' 2Tl and the later creation of 
a sea pair 5292 giving rise to the two mesons forming 
T 2 ■ {Qq2) miT2 ~ (Qq2) W2T2 ■ We shall call this two me- 
son threshold state |T 2 ) . 

Let us now imagine that we have two (or more) close 
thresholds in the sense that their effects overlap. In such 
a case we shall assume their net effect can be represented 
by one effective threhold whose properties will be ana- 
lyzed below. By using these effective thresholds we can 
generalize the form of V(r) by assuming that in between 
any two thresholds the same physical mechanism as de- 
scribed above between separate (non overlapping) thresh- 
olds is at work. 

For the sake of simplicity we shall reduce the size of 
the transition regions from one form of the potential to 



other just to the crossing points rr«- In order to unify 
the notation we define Mt = (note that To does not 
correspond to any real meson-meson threshold) so that 
the Energy Dependent Potential V(r) reads: 

1. M T „ < E spec < M Tl : V{r) = Vm Tq -m Ti M 
k 

Vm To -m Ti (r) = or - - r < r Tl 

Vm Tq -m Ti (r) = M Tl -m,Q- m-Q r > r Tl (1) 

2. M Ti < E spec < M Tz+1 : V(r) = Vu T . -M T . +1 (r) 
(* = 1,2...) 



V Mt . -M Ti+1 (r) = M Ti -m,Q- m-Q r < r Tz 
k 

Vm Ti -M Ti+1 (r) = or - - r Tt <r < r Tl+1 

V MT .-M Ti+1 (r) = M Ti+1 -mQ-m-Q r> r Ti+1 (2) 

with the crossign distances defined by 
k 

or Ti = M Ti - rriQ - m-Q- (3) 



This potential V(r) is drawn in Fig. [5]for the two lower 
energy regions. 

Let us now imagine that T n for example is an effective 
threshold resulting from the overlap of and T% that 
we shall call its parent thresholds. Then, in the over- 
lap region the probability for the system having a two 
static-light meson component has contributions from 
as well as from Otherwise said a given probability 
for having a two static-light meson component is reached 
at a shorter distance than in the case with only TjJ) or 
I T^j. Therefore we expect rj 1 ,, < fyi, r T 2 and conse- 
quently a mass for the effective threshold lower than the 
masses of its parent thresholds: Mt„ < M T i, M T 2. We 
shall comment about the physical interpretation of this 
prescription in the next section . 

Let us remark that the EQM interaction V(r) con- 
tains the sea quark pair effect through the coupling to 
thresholds. Therefore it corresponds to an unquenched 
quark-antiquark potential. It should be pointed out that 
the unquenching process implies in general the creation 
of any number of sea pairs and their recombination with 
the static (QQ) giving rise to many meson states. For 
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FIG. 2: Energy dependent potential V{r). The solid (dashed) 
line indicates the potential in the first (second) energy region 
for + (0 ++ ) bb states with m b = 4793 MeV (threshold masses 
and crossing radii from Table lU below) . 



instance the creation of two pairs q a qZ and qbgE can re- 
combine with (QQ) to generate two static-light mesons, 
(QOa) and (Qqb), plus a light-light meson (q a qb)- How- 
ever it may be reasonably thought that higher order pro- 
cesses like this involving the creation of more than one sea 
pair have less influence in the quark- antiquark interaction 
than the considered processes mediated by one sea pair 
creation. Otherwise said we have assumed that the un- 
quenched quark-antiquark interaction comes mainly de- 
termined by the two static-light meson thresholds. 

B. Threshold Effects 

The EQM interaction expresses the idea that the open- 
ing of a threshold introduces an energy dependence in the 
quark-antiquark potential. 

To be more precise let us imagine first that we had 
only one threshold T\ . Then the potential would be 

M To < E spec < M Tl : V(r) = V Mtq -m Ti (r) 

M Tl < E spec < M Too : V(r) = V Mti -m Too (r) (4) 

where Vm Tq -m Ti {?) is given by ([lj and we have defined 
Mt x = oo and Vm Ti —m Too ( r ) from (i = 1, 2...) 

Vm Ti -m Too (r) = M T% - vriQ - m-Q r < r Tz 
k 

V Mt . -m Too (r) = or r > r T% (5) 

r 

If we additionally define Vm Tq -m Txi {f) as the quenched 
(corresponding to the non threshold case) Cornell poten- 
tial (0 < r < oo) 

k 

Vm To -m Too (r) = V C or(r) = or (6) 



then it is clear from ((4]) that the effect of T\ is to intro- 
duce an energy dependence in Vcorif) = Vmt„-m Too (?*)■ 
It is also clear from the form of this energy dependence 
(the steep increase of the Cornell potential is cut at the 
threshold giving rise to Vm Tq -m Ti ( r )) that the opening of 
threshold T\ produces attraction (mass decrease) on the 
eigenstates of Vcor (r) with energy close to Mt x (the pre- 
cise numerical meaning of this attraction will be estab- 
lished later on), the closer in energy the more attraction 
on them. 

Let us now consider that we had only two thresholds 
T\ and Ti. Then the potential would be 

M Ta < E spec < M Tl : V(r) = V Mto -m Ti (r) 
M Tl < E spec < M T2 : V(r) = V Mti -m T2 {t) 
M T2 < E spec < M Too : V(r) = V Mt2 -m Tx (r) (7) 

By comparing Q and ([7]) it is clear that the opening of 
T 2 has no effect on the form of Vu T -M T ( r ) (the poten- 
tial below Mt x ) but introduces an energy dependence in 
Vm Ti -m Too (r) which causes attraction (mass decrease) on 
the eigenstates of Vm t -M Too ( r ) close in energy to Mt 2 - 

By substituting (Ti, T2) — > (Tj,Tj + i) we can generalize 
these results to the opening of any other threshold Tj+\ 
in a straightforward manner. 

From these considerations we may conclude that the 
eigenstates of Vm t ._ 1 -m t .( t ) are obtained from i) the 
eigenstates of Vm t ._ -m Tco with energy between Mx i _ 1 
and either not modified or attractively modified by 

the threshold Ti and ii) the eigenstates of Vm t . -m Too 
with energy above Mj< t attractively modified by the thresh- 
old Ti giving rise to states with energy below M^. 

To fix the ideas let us particularize this conclusion to 
the lower energy region: let us consider, for a definite 
set of quantum numbers I(J PC ), the quark-antiquark 
spectrum below where T\ is the first threshold for 
such set of quantum numbers. Then Vm To -M Too (r) = 
Vcor{r)- According to our analysis T\ causes attraction 
on the Cornell eigenstates with energy close below or 
above Mt x ■ As a consequence some Cornell eigenstates 
above threshold may be modified to acquire an energy 
below threshold. Therefore the EQM spectrum below T\ 
(E spec < Mt x ) is formed by the non modified Cornell 
eigenstates plus the attractively modified Cornell states 
lying below T\. 

From this discussion on threshold effects we may phys- 
ically interpret the prescription given in the former sec- 
tion for the mass of an effective threshold T n . Let be 
T,j and T% the two overlapping (parent) thresholds so 
that M T i < M T 2. Their joint effect makes the EQM 
corresponding eigenstates to have less energy than the 
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eigenstates of Vm t _ -M T 2( r ) since T* produces an ad- 
ditional attraction. This is equivalent to consider these 
EQM eigenstates as resulting from Vm t _ -« Tll (r) be- 
ing T„ an effective threshold with mass Mt„ < M T i . On 
the other hand the EQM corresponding eigenstates have 
less energy than the eigenstates of Vu Tl -M Tn+1 (?") since 
T% produces an additional attraction. This is equiva- 
lent to consider these EQM eigenstates as resulting from 
Vm t „ -m t „ x ( r ) being T n an effective threshold with mass 
Mx n < M T i . Putting together both prescriptions we get 
Mt„ < M T i , M T 2 as expected. 

C. Eigenstates 

Let us suppose a definite set of quantum numbers 
I(J PC ) for which we have a set of thresholds that we shall 
call simply {Tj} and a set of corresponding unquenched 
EQM eigenstates. From them let us choose an eigenstate 

of VM Tt l -M Tt (r) that we shall call (QQ)(fc) Ti _ 1 _ Ti ^ , 
where k stands for its spectroscopic identification (k = 
n r l in terms of the radial, n ri and orbital angular momen- 
tum, I, quantum numbers), with energy E^) T . T . G 

[M^.^MtJ . According to our discussion on threshold 
effects we may approximate this eigenstate, in first order 
perturbation theory, by 



(QQ) 



(k)T i _ 1 -T i 



{QQ) 



i k )T i _ 1 -T„ 



(8) 



where 



(QQ)(k) T - t-Too/ denotes an eigenstate of 
VMr-^-MT^i 1 ") below or close above the threshold Tj 
and |Tj) stands for the threshold state formed by two 
static-light mesons in a relative S— wave, |T) = 
\Tti Ti fBtxr t ) , with M Ti = M miT . + M m2T . . Henceforth 

we shall call 



(QQ)(k) T . _ Xoo / the generator state. 



As the generator 



and threshold IT,-) 



(QQ)(k) Tl _ 1 -: 

states are mutually orthogonal (note that in the i = 1 
case they correspond to Cornell and threshold eigen- 
states) the coefficients oi{k) T , T . and (3{k) T . T . satisfy 



( (t)T f _ r 



1 



(9) 



so that we can write them in terms of a mixing angle 

Q<0(*w 



< § in the form 



^(*)t 4 _ x -t ( = sine (k) Ti _ 1 - Ti (10) 

Let us take for instance the i = 1 case. Any EQM 
eigenstate (QQ)(k) To _ Tl ^> with energy less than Mt x cor- 
responds to an eigenstate of Vm t —m t ( r )- By compar- 
ing the forms of V MTo -M Tl {r) and Vm To -m Too (r) given 



respectively by (fT]) and © it is clear that if we project 
{QQ)(k) To -T 1 ^ to S e f a wave function in configuration 
space, its short range part (r <C ) will be determined 
by the generator state (QQ)(fe) To _ Too ^> ■ This argument 
can be immediately generalized to any other EQM eigen- 
state (note that on the contrary the r 3> behav- 
ior can not be correctly obtained from the approximate 
perturbative form of the eigenstate). Therefore the co- 
efficients «(fc) r . T . can be calculated by direct com- 
parison (division) of the short range wave functions of 



and 



(QQ)(k) Ti _ 1 - Ta 



Then the coef- 



(QQ)(k) Ti _ 1 - 

ficients P(k) T . _ T . can be obtained from the normaliza- 
tion condition ©. 

Hence we may conclude that the decomposition © 
tells us about the general spatial structure of a quark- 
antiquark EQM eigenstate: at short distances it corre- 
sponds to a conventional (resulting from a linear confin- 
ing model) Quark- Antiquark whereas at long distances it 
does to Meson-Meson. Obviously the Meson-Meson com- 
ponent has to correspond to a molecular state and not 
to the free threshold state we have used. Nonetheless 
we shall assume that the approximate decomposition (0) 
can be employed by treating |Tj) as a molecular state. We 
shall make use of this assumption for a simple qualitative 
analysis of the decays of an EQM eigenstate in terms of 
its generator and Threshold (molecular) states. 

Henceforth we shall refer to the EQM and Cornell 
eigenstates as EQM and Cornell states respectively. 



D. Strong Decay Modes 

Let us consider for a definite set of quantum numbers 
I(J PC ) a general EQM state [QQ){k) T ._ x _ T . \ with en- 
ergy -E(fc) T . T . above the first absolute threshold (note 

that this threshold may correspond to other J PC quan- 
tum numbers). 

Regarding strong decays we shall assume, in ac- 
cordance with our interaction model, that the domi- 
nant decay modes from the Quark- Antiguark component 

[QQ)(k) T . t _ T ^ correspond to transitions to open - 

flavor two meson final states. 

Actually this is an usual assumption when evaluating 
decays in conventional quark models although the quan- 
titative implementations of this mechanism carried out 
until now have not allowed for a precise computation of 
strong widths (see for instance d, Hj]). 

Let us consider first the i > 1 case. The dominant de- 
cay modes will be then the S— wave open flavor meson- 
meson channels with energy close below E^ T _ T . 
Then if T,*_i is a regular (non effective) threshold the 
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dominant mode will be the S— wave 97tiT i _ 1 97t2T;_i chan- 
nel. Else if Tj_i is an effective threshold, resulting from 
the two parent thresholds T}_ x and Tf_ x , and M T i , 
M T 2 < EnA T _ T , then the dominant decay modes 
will be the S— wave 97? 1 ' r i 1 9^2T. 1 1 an d 97T 1T 2 i 97t 2 T 2 
channels. Else if Tj_i is an effective threshold and 



E (k)T i ^ 1 -T z < M T}_ V 



Mrr 



then the dominant decay 



mode will be the 5— wave open flavor meson-meson chan- 
nel with mass closer below Mr 1 _ 1 (notice though that P 
or D~ waves could also contribute if there is not available 
S— wave channel or it is far below threshold; in this case 



even virtual S— wave 971 



it. 1 



971 



27? 



and 97t 



l'i',2 



971 



27? 



channels could be active if E^) T . Tt ~ Mr i _ 1 ). 



Let us examine now the i — 1 case. As P. 



< 



the dominant decay mode will be a P (or D) wave 
open flavor meson-meson channel with mass closer below 
Mt x - Of particular interest will be the case when T\ is 



a non effective threshold and E, 



(fc) To _ Tl 



M Tl . Then 



the dominant decay may proceed through the virtual S— 
wave 97tiT 1 97t2T 1 channel. 

As for the Threshold (molecular) component |Ti) it has 
been argued that due to the color structure of the quark- 
antiquark interaction (whose spectroscopic contribution 
is implicit in the values of the parameters used for the po- 
tential) the leading order meson-meson interaction, nec- 
essarily involving quark exchange between the mesons, 
results in a coupling with two body channels contain- 
ing a hidden flavor state {QQ) plus a light-light meson 
Since the energy of the EQM state, Ey.) T T ., is 
smaller than Mj,, it seems reasonable to consider that 
the \Ti) decay may proceed through those of these two 
body channels which are kinematically allowed. In par- 
ticular S— wave two body channels with mass close be- 
low £?(fe) T . T . can be expected to be favored (again P 
or D— waves could also contribute if there were not S— 
wave channel close below threshold; in this case even a 
virtual S— wave channel at about threshold could be ac- 
tive if B(fc) T . T . — MTi)- Moreover, as the Q mass is 
much bigger than the light quark one we may assign the 
orbital angular momentum to the light quarks so that 
the formation of an I = hidden flavor state may be 
expected to prevail. 

Therefore we shall assume for the |Tj) component that 
S— wave ((QQ)i=a+ light-light meson) channels coming 
out from quark exchange and with masses close below the 
energy of the EQM state are the dominant decay modes. 

As mentioned above we shall centre on the qualita- 
tive description of the decays. Nonetheless we shall try 
to combine results from conventional quark and molecu- 
lar models with our probabilistic description to extract 
from data rates (when available) some quantitative pre- 
dictions. 



E. El Transitions 

Cocerning El electromagnetic decays QQ — > 
l{QQ) final we shall assume Extended Vector Dom- 
inance (EVMD) 0. 

For the threshold component |Tj) EVMD implies the 
decay occurs dominantly via vector meson dominance 
(through p and tu conversion into a photon). 



For the (QQ)(k) T . 1 _ T ) component EVMD implies 

an Elementary Emission Model (EEM) such that the 
photon is emited by the Quark or the Antiquark. Actu- 
ally this is the usual mechanism assumed in the i = 1 case 
corresponding to the decay of Cornell states. In this case 
the application of nonrelativistic, zero recoil and dipole 
approximations makes clear the sensitivity of some tran- 
sitions to wave funcion details. 



Again we shall centre in a qualitative description of 
these processes although we shall try to generate some 
quantitative predictions from the application of the EQM 
to reproduce available data. 



F. Leptonic Decays 

QQ — > e + e~ transitions take place through a virtual 
photon. For 3 Si heavy quarkonia states to which we shall 
restrict our attention the transiton probability is deter- 
mined by the QQ wave function at the origin $gg(0). As 
the short range behavior of any EQM state corresponds 
to that of its generator the leptonic decay will proceed 
dominantly through the Quark- Antiquark component. 

Thus for a state below the first threshold, i = 1, the 
decay will proceed through its Cornell generating state. 
Unfortunately, due to uncontrolled QCD and relativistic 
corrections, no precise theoretical estimations of absolute 
leptonic widths for Cornell states have been achieved un- 
til now, only the ratios between leptonic widths approach 
the experimental values (see @ and references therein). 

Dimensionally the leptonic width depends on 



divided by a square mass factor. For a 



Cornell state this factor is usually approximated by the 
square mass of the state E? LA . For an i > 1 EQM 
state this factor may involve the mass of the eigenstate 



E. 



as well as the mass of the lower threshold 



-^7;_i that enters in the determination of the generator 
state. 

Taking these considerations into account we shall as- 
sume for the 



(QQ)(k) T . _ T . / leptonic width the general 
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[i > 1) form 



T[(QQ) {k)Tii _ 



*(QQ) 



{k) T 



(0) 



(11) 



M 2 (E, 



,M t , 



where the square mass factor M 2 must satisfy 
M 2 ( J B (fc)To _ Ti ,M To ) = ^ )ro _ ri . 

We shall make use of this assumption to try to fix 
the specific form of M 2 (E^ T . ,Mr j _ 1 ) from phe- 
nomenology. 

G. Decay Widths 

Let us consider a general decay in the form QQ — > AB. 



Let be 



(QQ) 



(fc)T 4 _ 1 -T j 



the initial state. The physi- 

(QQ)(k) Ti _ 1 - T , 



cal mechanisms giving AB from the 

and \Ti) components can be considered orthogonal in the 
sense that the application of the same decay mechanism 
to both components give rise to different final states. 
Think for example of the quark exchange process. It 
gives rise to |T,-) — > {QQ)i=o+ light-light meson. How- 



ever when applied to 



{QQ) 



it can not reach 



the same final state. Hence we shall write the width for 
a general decay as 



T([QQ) (k)Tt i _ T ^AB) ~ 



cos 2 (9 



12) 



the complexity deriving from overlapping thresholds has 
been parametrized in terms of artificial effective thresh- 
olds with undetermined masses. 

Anyhow, keeping in mind these possible shortcomings, 
we think it is worthwhile to examine the physical conse- 
quences deriving from this simple dynamical model for 
heavy quarkonia to try to learn from them possible av- 
enues for future progress. 

Aiming at a joint description of charmonium (cc) and 
bottomonium (bb) we shall use for both the same values 
for the parameters of the potential a and k. Let us realize 
that in the first spectral region E spec < Mt x , for energies 
not close to the threshold, we hardly expect any threshold 
effect. In other words Vcor{r) should describe reasonably 
well this part of the spectrum. Actually this is the case. 
It turns out that for a value of the Coulomb strength 
k = 100 MeV.fm corresponding to a strong quark-gluon 
coupling a s = || ~ 0.38 (in agreement with the value de- 
rived from QCD from the hyperfine splitting of lp states 
in bottomonium [Tl| and also with the value obtained 
from the fine structure splitting of lp states in charmo- 
nium [12]), one can choose correlated values of a and mq 
to get such description. In this regard, as we are deal- 
ing with a spin independent potential, we may compare 
as usual the calculated s— wave states with spin-triplets, 
the p— wave states with the centroid obtained from data 
and the d— wave states with the only existing experimen- 
tal candidates. Indeed it would be better a comparison 
with the centroids for all states but the dearth of spin 
singlet data makes this unfeasible. Henceforth we shall 
fix cr = 850 MeV/fm, m b = 4793 MeV and m c = 1348.5 
MeV so that the differences from the calculated Cornell 
masses to data below the first thresholds are less than 
30 MeV in bottomonium and 60 MeV in charmonium 
(these differences could be reduced when comparing with 
the centroids for all states). For m c its value has been 
fine tuned to get an EQM mass for X(3872) within its 
experimental energy range. 



H. Aplication to Heavy Quarkonia 

Certainly the EQM just defined might be too simplis- 
tic for spectroscopic purposes. First, the (quenched) Cor- 
nell potential used does not contain spin dependent terms 
that, apart from relativistic corrections, we know may be 
significant for the lower energy states. Second, only open 
flavor thresholds have been considered and no thresh- 
old widths have been taken into account. Third, the 
coupling between quark-antiquark and any threshold has 
been chosen with the same form when approximating it 
by an abrupt change in the potential at the crossing radii; 
however, it is possible that i) the creation of uu and dd 
pairs is favored against ss creation producing some dif- 
ference in the coupling to corresponding thresholds and 
ii) the coupling to a threshold with several possible val- 
ues of J is different for each of these values. Fourth, 



1. Thresholds 

In order to apply our model to a particular set of heavy 
quarkonia states with definite I(J P ) we have to look for 
meson {Qq) - meson {Qq) thresholds (q : u, d, s) coupling 
to these quantum numbers. We shall consider the two 
mesons to be in a relative S— wave so that the threshold 
energy corresponds to the sum of the masses of the con- 
stituent mesons. The lower thresholds for charmonium 
and bottomonium are listed in Tables U and UU 

It is important to remark that we have used isospin 
symmetry to construct thresholds with well defined 
isospin. This means that we are neglecting the mass 
differences between the electrically neutral and charged 
members of the same isospin multiplet, for example D 
and D ± with PDG quoted masses 1864.91 ± 0.17 
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and 1869.5 ± 0.4 or more dramatically Z?^(2400)° and 
Do (2400)- with quoted masses 2318 ± 29 MeV and 
2403 ± 14 ± 35 MeV respectively. Indeed we have cal- 
culated the threshold masses by taking for any member 
of the multiplet the lower mass value (1865 MeV and 
2318 MeV in the examples). We shall comment later on 
the consequences deriving from isospin breaking. 

Regarding C parity we rely on the quark model assign- 
ment of C parity values to charmed and bottom mesons 
(see for example the Quark Model Section in [13j) to ob- 
tain the C parity value for the threshold as the product 
of the C parities of the component mesons. Thus the sign 
preceeding the charge conjugate term (c.c.) is determined 
from the required threshold C parity. 

Unfortunately not all thresholds are experimentally 
well known. For example in Table U there is one thresh- 
old with D°D 1 (2430)°, which contains a 3 Pi - 1 P x 
mixing state Z?i(2430) with a quite uncertain mass 
2427 ± 26 ± 25 MeV. This threshold may overlap (or 
not) with the similar one containing the other state of 
the 1 P\ - 3 Pi mixing Di(2420). For the sake of sim- 
plicity (less number of separate thresholds) we have 
assumed they overlap and have substituted this pair 
of overlapping thresholds by an effective one with a 
mass smaller than the well established mass in the 
pair (M D0 D ( - 2420 - ) o)- More overlapping is present for 

thresholds L>*(2007)°L>i(2420)°, L>*(2007)°L>i(2430)°, 

£>+D s i(2460)- and £>*+D* (2317)~, with masses about 
4429 MeV. Again we shall consider all of them forming 
one effective threshold. 

A quite similar or even worst situation appears in bot- 
tomonium, see Table [TT1 where there is a known 1 P\ — 3 Pi 
mixing state Bi(5721)° but its 3 P 1 - 1 P 1 partner B 1 (?) 
is unknown. We expect this missing state to have a 
mass close to that of i?i (5721)° giving rise to an effective 
threshold. Besides there is a i?}(5732) with quantum 
numbers not established yet that we have tentatively as- 
signed to J p = + . 



2. Spherical well approximation 

We may have a quantitative idea of the effect produced 
by thresholds by realizing, from Tables Q] and HH that in 
the energy region in between any two of them, Ti and 
Tj-|-i with i > 1, the potential can be approximated by 
a spherical well since the ar term is very dominant and 
and TT i+1 differ at most by 0.25 fm. So we can expect 
the Vm t . -M Ti ± (r) potential to give similar results to the 
spherical well 

V Mt*-m Ti+1 0) = M T t -mQ-rriQ r < a M+i (13) 



Kfln) 



FIG. 3: Spherical well (dashed lines) representing the quark- 
antiquark interaction in between the first and second thresh- 
olds for 0+(0++) bb states with m b = 4793 MeV (threshold 
masses and crossing radii from Table HT)l . 



where 



TTi +rr i+1 



(14) 



Graphically this approximation looks like Fig. [3] 

We know the general condition to be satisfied for the 

existence of bound states in the spherical well potential. 

It reads: 



Pi 



Jl+1 (^Hf^) ■ 1+1 \ h J 

~ l PU m / ■ V 



(15) 



where 



Pi = V 2 M (Espec - M Ti ) 



(16) 



pn = J2fi (M Ti+1 - E. 



spec J 



with fi being the reduced mass of the quark- ant iquark 
system and E spec the energy, ji and hj 1 ^ stand for Spher- 
ical Bessel functions of the first and third kind respec- 
tively (see for example Q)- 

In the particular case of s— wave states the conditon 
([T5"j) reads 



tan 



Piai,i+i 



Pi 
Pn 



(17) 



and the spherical potential has n bound states if the fol- 
lowing inequality is satisfied 



< 



y/2n(MT i+1 - Ai^Ki+i (2n - 1)tt 



(18) 



For P and D— waves the presence of the centrifugal 
barrier reduces significantly the number of bound states. 
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Moreover the lowest s— wave bound state is above Mr 4 

h 2 2 

by an energy greater than — . Therefore there will 

not be any s— wave bound state in between two thresh- 
olds separated in mass by less than this amount of en- 
ergy. We may use this result to quantitatively define 
the meaning of close thresholds with overlapping effects 
on the quark-antiquark interaction introduced in Section 
II.A. 

We shall assume that two thresholds Ti and Tj+i over- 

t2 2 

lap when they are separated in mass by less than j^f- — ■ 

In such a case the thresholds Ti and Ti + i will not 
be considered as distinct but giving rise to one effec- 
tive threshold. For the sake of consistency it can be 
easily checked, from Table HI that the assumed overlap- 
ping thresholds involving Di(2420)° and £>i(2430)° re- 
spectively satisfy the required condition if the nominal 
masses for these charmed mesons are used. 



III. CHARMONIUM 

A list of the thresholds used in the calculation of the 
charmonium spectrum (the end of the list is imposed by 
the current partial or total lack of knowledge of further 
open flavor mesons) appears in Table IIIII where a sim- 
plified notation with respect to the one used in Table H] 
has been employed: a threshold has been denoted by 
the first meson-meson term entering in the 1 = linear 
combination which defines its content. Thus, the first 
0(0++) threshold in Table ^(D°lf + D+D~) is specified 
as D°D and so on. 

As mentioned before the masses of the effective thresh- 
olds are not determined. We have denoted these masses 
by a question mark followed by a colon and by our 
guessed value which have been chosen so as to reproduce 
the charmonium spectrum. Thus we have chosen two 
effective thresholds with masses 50 MeV smaller than 

M Z^ZM2420)° and M D.(2007)°IM2420) reS P eCtivel y- Let 

us note that the second effective threshold substitutes 
four overlapping thresholds whilst the first one substi- 
tutes only two. However it is plausible that the creation 
of uu and dd pairs is favored against that of ss so that the 
effect of the two thresholds containing charmed strange 
mesons is small as compared to that of the two thresholds 
containing charmed mesons. This may justify our choice 
of two effective thresholds with the same difference in 
mass with respect to the overlapping ones. In any case 
increasing (or decreasing) the mass difference for the sec- 
ond effective threshold by less than 30 MeV would not 
alter the resulting spectral pattern. Notice also that we 
could increase the effective mass of the first threshold 
by 30 MeV or even more and still reproduce the known 
charmonium spectrum but the mass for A (4260) would 
be badly overestimated. 



The charmonium masses from the EQM are listed in 
Tables QV] and [V] as M E qm and compared to data from 
13]. For the sake of completeness they are also com- 
pared to the results from the Cornell potential Vcor ■ We 
have made use of the spectroscopic notation introduced 
in Section II. C with some simplifications: in the lowest 
energy region we call the states n r l and for Cornell states 
above the first threshold we use (n r l). 

As can be checked a reasonable mass parallellism be- 
tween the EQM spectrum and cataloged neutral char- 
monium, including conventional and non conventional 
states, is obtained. The lower energy part of the spec- 
trum is described by the Cornell potential which ceases 
to be valid when approaching the first threshold. Above 
this threshold the spectrum is determined by the relative 
energy location of neighbour thresholds and the corre- 
sponding eigenstates always get a non-negligible thresh- 
old component probability which can be calculated from 
the generator states as explained in Section II. C. The 
mixing angles 6{k) T - _ T determining these probabilities 
in all energy regions are listed in Tables IVTl and IVTI1 where 
the energy intervals, the EQM masses and root mean 

1/2 

square (rms) radii, Meqm arid ( r2 ) E Q M , as well as the 
masses and rms radii of their generator states, Mqen 
and ( 7 ' 2 )q^j V 7 are also shown. 

A glance at Tables IVTl and [*VTT1 makes clear that mixing 
is in general very relevant for states above the first thresh- 
old (notice that = 45° corresponds to equal probabil- 
ity for the conventional Quark- Antiquark and Threshold 
components). This is also manifest from the compari- 
son of the masses and rms radii for the EQM and gen- 
erator states: the masses Meqm are smaller than the 
corresponding Mqen as a consequence of the attraction 

1 /2 

caused by threshold whereas the rms radii (f 2 ) E qm are 

1 /2 

bigger than (r 2 ^ GEN as a consequence of the non negligi- 
ble probability of the threshold component. The higher 
the mixing the larger these differences. 

In order to establish the strong decay modes of the 
EQM states we need to know, apart from their masses 
and related thresholds, the masses for the S— wave 
((cc); = o+ light-light meson) decay channels coupling to 
the same quantum numbers. These masses are tabulated 
in Table EUH 

We shall consider the resonances /, uj and 4> as mixed 
states containing u5, dd and ss components. There- 
fore thresholds formed from charmed mesons as well as 
from charmed strange mesons can decay into two body 
channels containing these resonances. In order to estab- 
lish the decay of the Threshold component of a EQM 
state we shall assume the dominance of the mode with 
mass closer below the corresponding EQM energy. Notice 
though that when this dominant mode involves promi- 
nently a light quark flavor not present in the Threshold 
component some alternative decay could be favored. For 
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instance the J/ip4> mode involves prominently ss while 
the threshold component of the (0 + (0 ++ ), 1p)t 4 -t 5 state 

given by D°£>(2550) involves uu and dd. So an alterna- 
tive decay to J/ipu> might be favored. 

By combining the contents of Tables I Villi and IIIII with 
those of Tables |VI] and |VlT] we obtain from the strong de- 
cay model (see Section II. D) the dominant modes which 
appear listed in Tables HXl and IXl Then the comparison 
of masses and decay modes with data makes feasible a 
sensible identification of EQM states with charmonium 
resonances. As a matter of fact this assignment has been 
implicitly used in Tables IIVI and [V] to locate any EQM 
state and its assigned resonance in the same row. 

It is worth to realize that relativistic corrections to 
the calculated masses are expected to be much more im- 
portant for states in the first energy region than for in- 
terthreshold energy states. The calculation of v/c for 
the 1" states gives a value around 0.5 in the first en- 
ergy region to be compared with a value around 0.12 in 
the interthreshold regions. This reduces the effectiveness 
of the parameters in these regions and perhaps improves 
the precision for the calculated masses with respect to 
those of the first energy region. 

Following the PDG notation [l3[ we shall name the 
_ (1 ) EQM charmonium states as iP{Meqm) and the 
+ (J ++ ) as Xcj(Meqm)- Let us analyze more in detail 
the resulting spectrum. 

A. First Energy Region : E sp ec < Mt x 

In the first energy region E spec < Mt x most states are 
very predominantly Cornell states for which we assume 
that the implementation of non considered relativistic 
and spin dependent corrections in the EQM may give, 
to a large extent, proper account of their masses and 
properties. In this regard we expect mass corrections to 
be more significant for the lower energy states. Mixing 
angles 9 in this region have values below 20° except for 
the more prominent non-predominant Cornell resonance, 
the 2p state with quantum numbers 1 ++ that we assign 
to A (3872) as explained below. 

More precisely for J PC = 1 the Is, 2s and Id 
states can be considered as (almost) pure Cornell states 
(0 ^$ 5°). For 3s and 2d with mixing angles 03 S = 12° 
and 8ss — 15° the masses and the absolute values of the 
wave function at the origin are slightly lower than for 
their generator Cornell states. This makes the leptonic 
widths for these EQM states to be quite the same as for 
their generators. Therefore the calculated ratios between 
leptonic widths for states below the first threshold agree 
reasonably well with data. Concerning strong decays and 
as a consequence of the mixing we expect for these EQM 
states an increase of the J / ipirir decay width with respect 
to their corresponding Cornell generators since this mode 



is dominant for the Threshold component (see Table ITXl) . 
Unfortunately we have only an upper bound for the width 
r(i/>(3s) -> J/V>7T7r) < 320 KeV (see Table 12 in @ and 
references therein) so that we can not extract any con- 
clusion about this expectation. 

Regarding J c = J ++ the lp states for J = 1,2 can 
be assigned to (almost) pure Cornell states whereas there 
is some mixing for J = 0. In this last case we expect the 
El decay width to "fJ/ip to be slightly smaller than for 
the generator (note that due to the low mass of the state 
the contribution from the Threshold component through 
vector meson dominance is quite suppressed). On the 
other hand the 2p state for J = 2 should show some 
effect from mixing, in particular as the EQM mass is 
above that of cuJ/ip we expect an increase of the decay 
width to LuJ/ip (with respect to that of the generator 
state) caused by the Threshold component. As for the 
2p state for J = 1 the mixing is most relevant as we 
analyze next. 

1. X(3872) 

In the EQM the Xcl (3871.6) = (1++, (2p) To _ Tl ), 
which we identify with A (3872), results from the at- 
traction produced by the D°D*(2007) threshold on the 
(1++, (2p) To -T oa ) Cornell state with mass 3911 MeV. As 
pointed out before the precise coincidence of the EQM 
mass with the experimental one has been required to fine 
tune m c in the hope that this will allow us to carry out 
a more precise probabilistic decomposition of this state 
in its Quark- Antiquark and Threshold components. As 
a matter of fact we could get any mass as close (below) 
as we wanted to the threshold by slightly changing m c . 
From Table EH] we have 6»L + t ~ 50° so that 

A(3872) - 0.64C c i(2p)+ 

+ 0.771- 1+ + (DD*) (19) 

where £ c and 1 stand for the Cornell and Threshold com- 
ponents respectively. 

From (ITOl) the decay rate to a given channel AB will 
be given by 

T(A(3872) ->• A) ~ 0.4r(C c i(2p) AB)+ 

+ 0.6r (* 1+ + (DW) AB) (20) 

Let us consider first the electromagnetic decays AB = 
jJ/ip and AB = Jip(2s). 

The width T (* 1+ + (DD*) -> Ji/j(2s)) is expected to 
be very small since (w, p) vector meson dominance is sup- 
pressed. On the other hand V (*i++ (DD*) -> jJ/ip) 
could be significantly smaller than T(( c i(2p) — > 7J/^0 
as it is the case for the molecular model considered in [J] 
(see Table 9 of this reference). If this were the case then 

r(A(3872) -> 7 y,(2s)) ^ r(C c i(2p) -» 7 V>(2s)) 
r(A(3872) -► 7J/VO r(C c i(2p) 1 ' 
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Unfortunately the left hand side ratio is not well mea- 
sured experimentally. Indeed measurements from differ- 
ent groups seem to be not compatible (see Table 12 in 
@ and references therein). By assuming that the cor- 
responding bottomonium state Xbifep) is almost a pure 
Cornell state (this will be proved later on) (xbi(2p) ~ 
Cbi(2p)) and recalling that the measured branching ra- 
tios to jJ/ip a- n d 7T(ls) are quite the same for Xci(lp) — 
Cci(lp) and X6i(lp) — Cbi(lp)) respectively we may ten- 
tatively assume 

r(Cd(2p) -> 7V>(2s)) _ r(Cw(2p) 7 T(2s)) 



From (J24J) and ([25]) we get 



r(C c i(2p) -> 7J/VO r(Cw(2p) -> 7 T(ls)) 



(22) 



where the right hand side is well measured experimentally 
131 



r(xbi(2p) 7 T(2s)) _ 0.199 ± 0.019 



= 2.2 ±0.4 (23) 



r(xbi(2p) -> 7 T(ls)) 0.092 ± 0.008 

Hence we may tentatively conclude that the A (3872) — >• 
jj/tp width is about half of the A(3872) -)• jip(2s) width 
although this result could be modified by the molecular 
contribution to the "fJ/ip decay. 

Let us consider now strong decays. As the Xfci(2p) 
wT(ls) decay has a branching ratio of 1.63% there may 
be also a comparable ratio for Cd(2p) —> us J ftp . From 
the estimated partial width for x&i(2p) — > wT(ls) : 
1.56 ± 0.46 KeV (see Table 36 in reference @) we may 
reasonably expect, using the parallelism between char- 
monium and bottomonium results, the partial width for 
Cci(2p) — > u J/ip to be of the same order of magnitude. 

We may also have a ^E' 1 ++ (DD*) — > uJ/ip decay 
due to the isospin breaking driven by the difference in 

mass between the D+L>*(2010)" and D°D*(2007)° com- 
ponents of the threshold. Actually the D+D*(2010)- 
mass (3880 MeV) coincides with M u + Mji^ mak- 
ing feasible the decay T (dTJ*) -)• 7r + 7r~7r° J/ip) 
through an intermediate uj. For the same rea- 
son (let us recall that D + D* (2010) ~ is a mixing of 
1 = and 1 = 1 states) it is also possible the 
^1++ (DD*) — > ir + ir~J/ip decay through an interme- 
diate p. The width r (DD*) -> tt+tt-tt J/ib) 
has been estimated with the molecular model of [4[ 
(see Table 8 of this reference) to be smaller than 
r (DlF) -> J/ip) . Although these quanti- 
tative estimations (720 KeV vs 1290 KeV) can not be 
taken for granted (see below) they seem to suggest 
that our estimated value for T (Cci(2f>) —tuij/ip) could 
be quite smaller than T (D~D*) -t 7r + 7r~7r° J/ip) . 
Therefore we may tentatively conclude that 



r(A(3872) -> tjj/tp) ~ 

o.er (dW) ->• wJ/i/S) 

A similar argument can be applied to Cci(2p) 
ir + Tr~J/ip so that 



(24) 



r(AT(3872) iT+ir-J/iP) 



0.6r (DD* 



tt + tt J/ip) 



(25) 



r(X(3872)->u>J/ip) 
r(X(3872)-+Tr + 7r- J/ip) ~ 

r(VK 1++ (£)^)^J/^) 
+ + (D'D 7 )^7T+n- J/ip) 



(26) 



The right hand side of this expression may be estimated 
from the molecular model of |4| (see Table 8 of this refer- 
ence) to be about 0.5 — 0.6 what may correspond to data 
within the experimental uncertainty (see Table 12 in Q). 

Furthermore we can use our previous results to estab- 
lish the approximate relation 



rpf(3872)-> 7 .//^) 
r(X(3872)->7r + 7r- J/ip) ~ 

0.4r(Cci(2p)->7J/^) 
0.6r(l' 1 + + (£>£>»)-Hr+7i— J/ip) 

Then we can use the rate (see Table 12 in Q) 

T(X(3872) -> jjjf) n?fi + ni 
T (A(3872) -> w+ir-J/i;) ~ 

to conclude that 

r(c c i(2 P )-^ 7 j/^)^ 

0.39r(* 1+ + (DD*) -> TT+n-J/iP) 



(27) 



(28) 



(29) 



although this result could be modified by the molecular 
contribution to the jJ/ip decay. 

Finally let us examine the A (3872) -> L>°D*(2007)° 
process. From data this decay is an order of magni- 
tude more prevalent than A (3872) — > n + ir~J/'ip (Q 
and references therein). According to our model the 

decay to D°D* (2007) should be the dominant mode 
for C lC x(2p). Then if we assumed the qualitative re- 
sult from the molecular model of reference |H that the 

(DD*) -> D°D*(2007)° width is quite smaller than 
the ^>i++ (DD*) — > it + it~ J/ip width we might conclude 

that the A" (3872) -)■ D°D* (2007)° decay proceeds mainly 

through Cci(2p) -)■ D°D*(2007)°. Taking into account 
the experimental value for the total width of the A (3872), 
Ftotai — 1-3 ± 0.6 MeV, this would imply 

T (Cd(2p) -> L>°L>*(2007)°) < 4.7 MeV (30) 

what is far below the quantitative result from some model 
calculations @, Q but it might be compatible with the 
Cornell coupled-channel model of reference [l5[ once the 
mass for the Cci(2p) state is fixed at 3871.6 MeV (see also 
Table 4 in Q for a comparative study of several models) . 

Putting together all the results obtained we may guess 
some quantitative intervals of values for the partial 
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widths: 

rpf(3872) -»• jJ/i/j) : 10 - 40 KeV 
T(X(3872) -> 7<0(2s)) : 20 - 80 KeV 

T (V(3872) -> vr+TT- J/0) : 40 - 160 KeV 
T (V(3872) -> wJ/0) : 20 - 80 KeV 

T (x(3872) -> L> D*(2007)°) < 1.5 MeV 

It should be remarked that the guessed intervals 
for electromagnetic transitions are quite compatible 
with the values calculated in reference [HI ] whereas 
the open flavor strong decay width may be ob- 
tained from the Cornell coupled-channel model of ref- 
erence [HI as explained above. Moreover the rate 
r (Jf (3872) -> ojJ/4>) /r (V(3872) -> tt+tt" J/0) is con- 
sistent with the one coming out from the molecular model 
of reference [|| . 

B. Inter Threshold Resonances (ITR) 

In the energy regions between two neighbor thresholds 
Tj and Tj+i with i > 1, most resonances are generated 
from the attraction caused by the upper threshold Tj+i 
on the lower energy eigenstate of Vm t . —Mt ( r ) which 
can be below or above TJj+i. Sometimes, as it is the case 
in the X2 — T3 region for ++ and 1 ++ , this eigenstate is 
far above the upper threshold (T3) so that the attraction 
is not sufficient to generate any (\p) T ^_ T ^ ITR (see Ta- 
bles |V] and |VII|- The same mechanism explains the shock- 
ing absence of d— wave 1 resonances from 4-200 MeV 
to 4100 MeV. Despite this and contrary to conventional 
quark models the EQM predicts more ITR than exper- 
imentally extracted resonances. In particular three 1~~ 
missing resonances that we denote as '0(4496), ■0(4537) 
and 0(4604) and two 1 ++ missing resonances that we 
denote as x c i(4003) and x c i(4190) are predicted. The 
experimental extraction of these resonances would be a 
definite proof of the validity of the EQM model. 

With respect to the PDG cataloged neutral X reso- 
nances apart from V(3872) all of them can be reasonably 
assigned to calculated 1 and J ++ ITR's with the ex- 
ception of V(3940). Let us examine in more detail these 
assignments to identify possible additional missing states. 

1. X(3915) 

The mass (3915.5 MeV) and the observed decay mode 
J/ifiuj of the V(3915) suggest its identification with 
the (0+(0++), lp) Tl -T 2 ITR at 3898 MeV. Alternatively 
V(3915) could be the same resonance as Xcii^p) (mass 
3929 MeV) to which we have assigned the (0 + (2++), 2p) 
state at 3903 MeV in the first energy region. If this al- 
ternative were correct then the (0 + (0 ++ ), 1p)t 1 -t 2 ITR 
would correspond to a missing resonance x c o(3898). 
Therefore it is plausible the identification of V(3915) 
with either a Xco or a \c2 state. 



2. X(3940) 

The decay mode for this PDG resonance (DD* seen, 
DD not seen, J/tpuj not seen) does not fit with any of 
the calculated ITR. The fact that X(3940) is produced 
in double charm production as well as f] c (ls) and ?y c (2s) 
suggests its identification with a r\ c resonance with quan- 
tum numbers + (0 From conventional quark mod- 
els the mass of ?7 c (3s) is expected to be about 40 MeV 
smaller than the mass of 7/; (4040) [| ; moreover addi- 
tional relativistic corrections might push down this value 
In the EQM the first S- wave 0(0"+) threshold 

for n c , D°Z^(2400)°, at 4183 MeV is an additional (al- 
though weak) source of attraction. Besides there is a 
0(1+") DD* threshold at 3872 MeV (differing from that 
111 Table U in the sign in front of the complex conjugate 
component (c.c.)) so that the dominant EQM decay mode 
of this n c is (DD*)p whereas DD and J/i/jlj can not 
be reached by quantum numbers. Putting all these re- 
sults together it seems quite plausible the identification 
of X (3940) with 7i c (3940). 



3. X(4U0), X(4160) 

The PDG mass (4143 MeV) and decay mode 
(J/t/xf)) of V(4140) suggest its identification with ei- 
ther Xco(4140) = (0+(0++),lp) T3 -T 4 or x c2 (4140) = 
(0 + (2 ++ ), lp)T!-To ■ Exactly the same situation occurs 
for X(4160) (mass 4156 MeV). Therefore the two ITR, 
Xco and \c2 can be assigned to these two resonances al- 
though their precise one to one identification is not pos- 
sible yet since it requires the experimental determination 
of the total angular momentum of at least one of the two 
resonances. 



4. X(4350) 

The lack of predictions for + (l,2 ++ ) states above 
their last well known threshold at 4224 MeV does 
not allow for a clear identification of this resonance 
(mass 4350.6 MeV) in spite of the fact that there is 
a (0+(0++), lp) T4 -T 5 ITR at 4325 MeV with a domi- 
nant decay to J /ip<j) as experimentally observed. Indeed 
the presence of this ITR suggests that 0+(l, 2++) ITR's 
from still unknown thresholds may also be close in en- 
ergy. Hence we prefer to leave open the assignment until 
more data on thresholds become available. 



5. Xd (4003) and Xd (4190) 

The analysis carried out until now allows for a rea- 
sonable assignment of an EQM state to any PDG listed 
neutral J++ resonance up to 4350 MeV. But the EQM 
predicts the existence of two more ITR in this energy 
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region, the Xei(4003)^ (0+(l++), lp) Tl -T 2 with pre- 
dicted decay modes, DD* and J/tpcj and the Xci(4190) = 
(0+(l++), lp) T3 _T 4 decaying to D S D^ and J/ip(f). We as- 
sign them to missing resonances. 



6. X(4260) and (4415) 

In the EQM the existence of the 0(4280) = 
(0~(1 ),1s)ti-t 2 ITR below the mass of first well 
known threshold L>L>i(2420) (4287 MeV) relies on the 
presence of an effective threshold at 42 37 MeV r esult- 
ing from the overlap of DL>i(2420) and DL>i(2430). This 
ITR may be identified with V(4260) (actually the choice 
of a lower mass for the effective threshold could fit pre- 
cisely the experimental mass 4263 MeV) since the seen 
decays J/0(7T7r, KK) correspond to dominant ITR modes 
through J/-0/o(98O). Notice that these modes come out 
from the regular Threshold component. Concerning the 
predicted modes (D* S D* S , D S D* S , D*TF) P from the Quark- 
Antiquark component the experimental siuation is much 
more uncertain what may be indicating the physical com- 
plexity associated to the lower effective threshold. 

The other ITR involving an effective threshold is the 
(0+(l — ),1s)t 3 -t 4 state at 4418 MeV which may be 
identified with -0(4415) (mass 4421 MeV) if we take 
into acco unt that t he S— wave dominant mode in Ta- 
ble |IX] D*Dl (2400) is quite below the EQM mass what 
may do the D— wave decay channel (Z?£>2(2460))£> to 
be dominant as experimentally observed. Regarding the 
Threshold component the expected dominant decays are 

tO J/0(7T7T, 47r). 

The 0(4415) — > e + e~ width has also been measured 
(0.58 ± 0.07 KeV). If we apply, following Section II.F, 
the usual approximate expression for the width (see for 
instance Q and references therein) modified by a square 
mass denominator M 2 (E^i s - jT3 _ Ti , Mt 3 ) ■ 

r^(4415)^e+e- — 4e^O! 2 (l - -gf*). 

|2 (31) 

J^(1 3 )T3-T 4 (°)| 

■M2(£ (ls)T3 _ T4 ,M T3 ) 

where eQ is the quark charge e c = 2/3, a ~ the 
electromagnetic coupling constant, a s the strong cou- 
pling strength and R(is) T3 _ Ti (0) the EQM radial wave 
function at the origin, then we can get a close value 
(0.69 KeV) to the experimental number by choosing 
M 2 (E {1s)t3 _ T4 ,Mt 3 ) = Ef ls)T3 _ Ti - Mf 3 . Unfortunately 
we have not more leptonic width data for interthreshold 
resonances to check the possible validity of this prescrip- 
tion. 



7. X(4360) 

The calculated mass and decay modes (D*Dw, 
0(2s)7T7r) for the ITR 0(4360) = (0-(l~ ), ls) T3 _ T3 
strongly suggest its assignment to V(4360) (mass 4361 
MeV and observed decays to D°D*~ir + and 0(2s)7t + 7t~). 

Notice that for the lower threshold T 2 = £>*£>* (2400) 
there may be a significant isospin breaking since 
D*(2400)° and D*(2400) ± with quoted masses 2318 ±29 
MeV and 2403 ± 14 ± 35 MeV respectively may sig- 
nificantly differ in mass. Actually the observed decay 
D*-D°n+ points out to a decay through D*"L»5(2400)+ 
what could be indicating that the mass of this channel 
is pretty close below the mass of the ITR. This would 
mean that the mass of the Dq(2400) ± would be close to 
the lower value of its experimental mass interval (2354 
MeV). 



8. ip (4496), i> (4537) and i> (4604) 

Three missing resonances, 0(4496), 0(4537) and 
0(4604) are predicted in between the experimentally 
known 0(4415) and V(4660). They correspond to the 
(0-(l~),ls) Ti _ Ti+1 ITR's with i = 4,5,6. For 0(4604) 
the predicted mass and dominant decays may be altered 

by the presence of two still unknown thresholds D°Di(?) 
since two D\ mixing states ( 3 Pi — 1 Pi) with masses 
around 2800 MeV can be expected (there is a Z?(2750) 
with non established quantum numbers that could well 
be a + sate). For 0(4496) and 0(4537) their vicinity 
may make difficult their experimental disentanglement. 
In this respect their differentiated dominant decays could 
help to their extraction. 



9. X(4660) 

The V(4660) (mass 4664 MeV) with seen decay to 
0(2s)7r+7T- is identified with the _ (1 — ), 1s) t ?-t s ITR 
at 4673 MeV although the open flavor decay of V(4660) 
to D°D*~ir + docs not correspond in principle to the 
predicted dominant decay (D*D*K) of the ITR. How- 
ever, as explained for 0(4604) there may be an effec- 
tive threshold from D°Di(?) with mass close to that of 
£>*+£> s i(2536)-. If this occurred the ITR mass could 
still be within the experimental mass interval but with 
the dominant open flavor decay to DD*ir as observed. 



C. Bottomonium 

The list of thresholds employed in the calculation of the 
bottomonium spectrum appears in Table IXI1 The lack of 
knowledge about further thresholds prevents extending 
the list to higher energies. 
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As it was the case for charmonium there may be two 
possible effective thresholds. They may correspond to the 

pairs of close thresholds (B° J Bi(5721)°, B°B 1 (?)°) and 

( J B*°Si(5721)°, B*°Bi(?)°). From our study of the spher- 
ical well approximation it is clear that the overlapping 
energy interval between two thresholds is shorter in bot- 
tomonium than in charmonium since it depends on the 
inverse of the reduced mass of the system. Therefore we 
expect the mass for the effective thresholds in bottomo- 
nium to be closer to the mass of the parent thresholds 
as compared to charmonium. So we have chosen two 
effective thresholds with masses 25 MeV smaller than 
M b° b i( 5721) ° and M b«° b 1 (5721) ° respectively since they 
may reproduce data as we shall see. As a result the 
upper effective threshold is located at about the same 
energy that the first regular threshold (11023 MeV). We 
should note though that this regular threshold does not 
overlap with the parent thresholds of the effective one. 
In practice this means that we shall only consider one 
threshold at 11023 MeV. 

The bottomonium masses from the EQM are listed in 
Tables IXlTl and IXIIII and compared to data from [l3| and 
to results from the Cornell potential (notation as in Ta- 
bles |T3 and E| . 

At difference with charmonium unconventional reso- 
nances have not been clearly identified experimentally in 
bottomonium until now. According to the EQM the only 
PDG cataloged neutral resonance that may correspond 
to an ITR is T(11020). For the rest of PDG resonances 
an unambiguous mass parallelism with calculated EQM 
states can be established. 

The generator masses and mixing angles for the EQM 
states appear in Tables IXIVI and fXVI (notation as in Ta- 
bles IVII and IVII[) . The same comments done in charmo- 
nium about the relevance of mixing above the first thresh- 
old can be traslated to bottomonium. From the mixing 
angles the Quark- Antiquark and Threshold components 
probabilities can be calculated. 

In order to assign the EQM states to PDG resonances 
we proceed to an identification of masses and dominant 
decay modes. For this purpose the S— wave ((bb)i = o+ 
light-light meson) decay channels are listed in Table IXVfl 
and the dominant strong decay modes from the Quark- 
Antiquark and Threshold components of the EQM states 
are compiled in Tables IXVIII and IXVIIII As in the char- 
monium case we shall consider the resonances /, u and 
(f> as mixed states containing iiu, dd and ss components 
in order to establish the dominant decay mode of a EQM 
state from the channels with mass closer below the cor- 
responding EQM energy. Notice though that when the 
dominant mode involves prominently a light quark flavor 
not present in the Threshold component some alternative 
decay could be favored. For instance the T(ls)0 mode 



involves prominently an ss while the Threshold compo- 
nent of the (0+(0_++), lp) Tl -T 2 state given by B*°B~*° 
involves uu and dd. So we have also listed an alternative 
decay to T(ls)w. 

As it happened in charmonium the calculated values 
of v/c for the 1 states are much bigger in the first en- 
ergy region (0.3) than in the interthreshold ones (0.05) 
where the effect of relativistic corrections get reduced. 
Therefore we may expect the calculated masses in the 
interthreshold regions to be more precise than the calcu- 
lated ones in the first energy region. 

Following the PDG notation we shall name the 
_ (1 __ ) EQM bottomonium states as T(Meqm) and 
the + (J ++ ) as xw(^£Qm)- Let us analyze more in de- 
tail the resulting spectrum. 



D. First Energy Region : E sp ec < Mt-l 

In the first energy region E spec < Mj^ most states are 
almost pure Cornell states (8 < 5°) for which we assume 
that the implementation of non considered corrections in 
the EQM may give, to a large extent, proper account of 
their masses and properties. A signature of the small 
mixing for the J ++ 2p states (8:2 — 3°) may be the 
surprisingly large Xfci,2(2p) — > wT(ls) transition proba- 
bility experimentally observed (see Table 36 in Q and 
references therein). Despite the small Threshold compo- 
nent probability (4%) the fact that the wT(ls) is the only 
dominant channel kinematically available for this compo- 
nent being the wT(ls) mass very close below the EQM 
mass might explain the observed rate. 

The J ++ 3p states present a significant mixing, par- 
ticularly the J — state (8 3p (0 ++ ) = 24°). Therefore we 
expect for these EQM states an increase of the wT(ls) 
decay width with respect to their corresponding Cornell 
generators. Regarding El decays the mixing effect de- 
pends on the contribution from the Threshold component 
through vector meson dominance. As this contribution 
is suppressed for Xbj(3p) —> 7T(2s) we expect a decrease 
of the 7T(2s) width with respect to that of their corre- 
sponding Cornell generators. 

For J PC = 1~~ only the 5s and Ad states (9^ s ~ 
10°, = 14°) have a significant mixing in the first en- 
ergy region. This has hardly any effect on the leptonic 
widths whose values are quite close to the Cornell ones 
(the EQM masses and the absolute values of the wave 
function at the origin are slightly lower than for their gen- 
erator Cornell states what produces no variation in the 
calculation of the leptonic widths) . Hence the calculated 
ratios between leptonic widths in the first energy region 
agree reasonably well with data. Mixing is however rel- 
evant for strong decays. A look at Tables IXlTl and IXVIII 
tells us that in particular the 5s and Ad decay widths to 
T(ls)(7T7r) and T(3s)(7T7t) can be significantly enhanced 
with respect to the Cornell case since the dominant two 
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body decay channels T(ls)/o,i,2(^ 1300,^ 1400,) and 
T(3s)/o(~ 500) have masses pretty close to the EQM 
masses. This may explain the huge signals observed for 
these decays around 10870 MeV (see Table 36 in (gj and 
references therein). Regarding the even higher signal to 
T(2s)(7T7r) we shall comment below. 

E. Inter Threshold Resonances (ITR) 

As it occurred in charmonium sometimes the attrac- 
tion caused on the generator state by the upper thresh- 
old is not sufficient to generate an ITR. This explains the 
absence of a 1 ++ ITR in between the lower thresholds T\ 
and T2 in Table IXim as well as the absence of d— wave 
1— ITR's bet ween % and T l+1 for i = 1,2 in Table ETj] 
(see also Tables 1x171 and KYj . 

For 1" the EQM predicts five ITR states between 
10995 MeV and 11195 MeV. As the first two, at 10995 
MeV and 11039 MeV, might be hidden within T(11020) 
(see below) the extraction of 1 resonances around 
11090 MeV, 11130 MeV and 11140 MeV would give 
strong support to the EQM although their proximity may 
cause some overlap making diffcult their experimental 
disentanglement. Actually this energy region has been 
experimentally scanned but a thorough analysis of data 
may be required. Additionally the EQM predicts a ++ 
ITR around 10620 MeV and 0, 1, 2++ ITR's around 10707 
MeV and 10813 MeV which should also be found for fur- 
ther support. 

Let us analyze in more detail these predicted ITR and 
their comparison to currently existing PDG resonances. 

1. T(11020) 

The quoted mass of this PDG resonance (11019 ± 8) 
MeV lies in between the (0~(1 — ),1s) Ti -t 2 ITR with 
a mass of 10995 MeV and the (0~(1 — ), 1s)t 2 -t 3 ITR 
with a mass of 11039 MeV. Although these ITR masses 
are related to a somehow arbitrary choice of the masses 
for the two effective thresholds it is for sure in the EQM 
the presence of at least two ITR, one below the well es- 
tablished threshold at 11023 MeV and other between this 
threshold and the well established one at 11072 MeV. 
This suggests that T(11020) could actually come out 
from the overlapping of these two ITR T(10995) and 
T(11039). Indeed the BaBar collaboration has reported a 
resonance with extracted mass 10996 ±2 MeV 18] which 
might be assigned to T(10995). Moreover a dominant 
decay mode from the Threshold component of this ITR 
is (T(2s)/o(980))„i rtua i (notice that the /o(980) mass un- 
certainty, 990 ± 20 MeV allows for this decay) what may 
contribute to explain the huge T(2s)7r7r signal observed 
around 10870 MeV. In this regard the T(11039) could be 
also contributing. Therefore we may conclude that there 
are several indications pointing out the validity of these 
EQM predictions although no definite conclusion should 



be extracted until the experimental confirmation or refu- 
tation of the existence of T(11039) as a differentiated 
resonance from T(10995). 

2. T(11090), T(11130) and T(11140) 

The EQM predicts three 1 — - missing resonances be- 
tween 10995 MeV and 11195 MeV. These are T(11090) = 
(0-(l — ), 1s)t 3 -t 4 , T(11130) = (0-(l — ), ld) T3 -T 4 and 
T(11140) = (0-(1~),2s)t 3 -T4- The T(11090) may be 
the best candidate for its experimental observation. For 
it we expect dominant decays to B* B*tt, B* Bit and to 
T(4s)7T7r and T{1s)(ljlj, 4tt, KK) what may help to its 
extraction. 



3. xw(10620) 

In the J++ sector the EQM predicts a Xbo (10620) = 
(0 + (0 ++ ), 1p)ti-t 2 we ll isolated resonance what may fa- 
cilitate its experimental extraction through the dominant 
decay modes BB and T(ls)w. 

4- Xi-0,1,2 (10704, 10708, 10710) and 
Xbo,i,2 (10784, 10815, 10822) 

Two mass triplets of XbJ resonances around 1707 MeV 
and 1813 MeV respectively are predicted. This has to do 
with the presence of common thresholds in the 0, 1, 2 ++ 
channels. Nonetheless both triplets are different in the 
sense that in the first triplet the three states share the 
same dominant decay mode from the Quark-Antiquark 
component (B*B*) whereas in the second one they do not 
being the dominant decay mode different for each mem- 
ber of the triplet {B S ~B~ S , B*B~f and B*W respectively) 
what may help to their experimental disentanglement. 



IV. SUMMARY 

From an educated guess for the ground state energy 
of a static quark - static antiquark system based on lat- 
tice calculations we have proposed an energy dependent 
quark-antiquark potential for calculating the spectra of 
heavy quarkonia within a nonrelativistic quark model 
framework. This potential, incorporating the interaction 
with open flavor meson-meson channels, allows for an 
unquenching of the conventional quark model that we 
call Extended Quark Model (EQM). Instead of the con- 
tinuous increase of the (confining) potential character- 
izing conventional quark models the energy dependent 
potential tends to a constant in any finite energy region. 
The EQM quark-antiquark eigenstates can be interpreted 
in terms of conventional (quark-antiquark) + (meson- 
meson) states. The conventional (quark-antiquark) com- 
ponent determines the short range spatial configuration 
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whilst the (meson-meson) component dominates the very 
long range part. The resulting physical image, in terms of 
conventional states, is that of an extended core of quarks 
surrounded by a mesonic cloud. This image has been 
used to get a qualitative understanding of the EQM decay 
processes from the dominant decay modes of the compo- 
nents. 

The EQM has been applied to the description of char- 
monium and bottomonium. The use of the same poten- 
tial in both cases is justified once the arbitrary constant 
in the potential has been absorbed in the definition of the 
quark mass. Regarding the level of precision of the cal- 
culated masses we may distinguish the first energy region 
from the upper ones since the form of the potential makes 
the inter threshold resonances (ITR) to become truly non 
relativistic states. Therefore we may expect an improve- 
ment of the precision for the calculated ITR masses as 
compared to the masses calculated for the almost Cornell 
states in the first energy region. More specifically each 
PDG cataloged cc and bb neutral meson has been reason- 
ably assigned to an EQM state with calculated mass close 
to the experimental one and with dominant decay modes 
fully compatible with the observed decay channels. Be- 
sides the EQM predicts the existence of non cataloged 
cc and bb neutral resonances whose discovery would give 
definite support to the model. In charmonium the pre- 
dicted absence of d— wave 1 resonances from 4200 
McV to 4700 MeV should also be emphasized. 

Certainly the model relies on threshold states whose 
meson components should be also consistently described 
by the EQM. In this respect the minimal form of the Cor- 
nell potential we have employed could not be sufficient 
for an approximate description of these mesons contain- 



ing light quarks and some refinements might be needed. 
Anyway there are some clear indications that a refined 
EQM should be considered for an appropriate descrip- 
tion of all meson sectors. For instance in the charmed 
strange sector the mesons D s o(2317) and D s i(2460) have 
unexpected low masses as compared to most conventional 
quark model calculations. In the EQM £) s o(2317) is a 
I(J P ) = 0(0+) state with mass below the first 0(0+) 
threshold DK {M DK = 2360 MeV) whereas L> s i(2460) 
is a 0(1 + ) state with mass below the first 0(1 + ) thresh- 
old D*K {M D , K = 2502 MeV). It is then plausible that 
these two states result from the attraction caused by the 
respective thresholds on the corresponding conventional 
(generator) states. In other meson sectors like the light 
unflavored one the EQM may also play a very relevant 
role in the understanding of the puzzle concerning 3 Pq 
states such as /o(980) lying close below the KK thresh- 
old. Furthermore the EQM generalization to baryons 
may provide a general scheme for the solution of endemic 
problems related to the description of some baryonic res- 
onances as AT* (1440), A(1405), A 5/2 _(1930)... 

In conclusion we have developed an unquenched quark 
model that may allow for a general description of hadrons 
as made of constituent quarks and antiquarks {qq for a 
meson and 3q for a baryon). The results obtained for 
heavy quarkonia with the simplest version of this model 
seem to point out that it incorporates essential physical 
ingredients needed for an accurate study of the hadronic 
structure. 
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Charmonium Thresholds Mt, VT t 
T\ (MeV) (fm) 



0(0++) 

D 0d° + D+D- 
DfD~ 

0(1 ++ ) 

{D°D* (2007)°+ 
L> + /J*(2010)~) - c.c. 



(D*(2007)°Z>(2007) + 
_D*(2010) + D*(2010)~) 

0(2++) 

(D* (2007)°D*(2007)°+ 
_D*(2010)+D*(2010)") 

Dt+DT 

0(1") 

(_D°Di(2420)° + 
D+Di (2420)") + c.c. 



(D°D 1 (2430) + 
D+D 1 ~(2430)) - c.c. 



(D*(2007)°L>5(2400) + 
D* (2010) + /J>5 (2400)" ) + c.c. 



(D*(2007)°L>i(2420) + 
L>*(2010) + Di(2420)") - c.c. 



(D*(2007)°L>i(2430) + 
D*(2010) + Di(2430)~) + c.c. < 4429 < 2.09 

D^D s i(2460)- - c.c. 

D; + D* 30 (2317)- +c.c. 

TABLE I: Meson-meson thresholds for cc states. Threshold 
masses calculated from the charmed and charmed strange me- 
son masses quoted in pl| . Inequalities have been used for the 
masses and crossing distances of the possible effective thresh- 
olds. 



3730 1.31 

3937 1.54 

3872 1.46 

4014 1.62 

4014 1.62 

4224 1.86 

< 4287 < 1.93 

4325 1.98 



ri tPC\ Bottomonium Thresholds Mt, 

Ti (MeV) (fm) 



0(0 +H 



0(1 



0(2"' 



0(1~) 



B°lf + B + B~ 10558 1.24 

B*°W° + B* + B*~ 10650 1.34 



(B°B* +B + B* )-c.c. 10604 1.29 
B*°W° + B* + B*~ 10650 1.34 



B*°B*° + B* + B*~ 10650 1.34 
BIB* 10830 1.54 



(B°Bi(5721) + 
S+Bi (5721)") +c.c. 



< 11003 < 1.73 



(B°Bi(?) + 
B+Bi (?)-)-& 



(£T° J B *(5732) + 
B* + BJ(5732)") + c.c. 



11023 1.76 



(B*°Bi(5721) + 
S*+Bi (5271) - ) - ex. 



< 11049 < 1.79 



[B*°Bi (?) + 
B* + Bi(?)-) + c.c. 



(B*°£ 2 *(5747) + 
B*+BJ (5747)") +c.c. 



11072 1.81 



TABLE II: Meson-meson thresholds for 66 states. Threshold 
masses calculated from the bottom and bottom strange meson 
masses quoted in For B}(5732) we have assumed J = 0. 
A question mark has been used for the mass of the unknown 
meson and inequalities for the masses and crossing distances 
of possible effective thresholds. 
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Ti 


Mesonl — Meson2 








21 




(o- 




? : 4237 




D°Di(2430)° 




T 2 


D*(2007) D5(2400) 


(1- 




4325 




i^zuu i ) u\ { ZQzyj ) 


I 1 






± 3 


D*(2007)°Di(2430)° 


(1- 




? : 4379 




D+D s i(2460)- 


(o- 






D*+D* 30 (2317)- 


(1- 


,o+) 




T 4 


D* (2007)° (2460) ° 


(1- 




4470 


r 5 


D+D s i(2536)- 


(o- 




4504 


r 6 


D s *+D s i(2460)- 


(1- 




4572 


r 7 


D*+D s i(2536)- 


(1- 




4648 


T 8 


D s *+D s * 2 (2573)- 


(1- 


, 2+ ) 


4685 



T 


D°D° 


(o- 


■-0-) 


3730 


T 2 


DfD~ 


(o- 


-.0") 


3937 


T 3 


TT(2007)°TT(2007)° 


(1- 


M-) 


4014 


Ti 


D*+D* s ~ 


(1- 


M-) 


4224 


T 5 


D D(2550)° 


(o- 


-.0") 


4405 



1++ 





D D*(2007)° 
TT(2007)°Z>(2007)° 


(o- 


M") 


3872 


T 2 


(1- 


M-) 


4014 


T 3 


D+D*- 


(o- 


M") 


4080 


T 4 


D*+D*~ 


(1- 


M-) 


4224 



Ti D*(2007)°D*(2007) (l~,l~) 4014 
T 2 D*+D*~ i 1 ', 1 ) 4224 



TABLE III: Meson-meson thresholds for / = cc states. 



J PC Notation MeQ 







MeV 


MeV 


MeV 


Is 




3046 


3096.916 ±0.011 


3046 


2s 




3632 


3686.09 ± 0.04 


3632 


Id 




3743 


3775.2 ± 1.7 


3743 


3s 




4061 


4039 ± 1 


4066 


2d 




4137 


4153 ± 3 


4143 




■r 2 


4280 


4263tg 




(1s)t 2 - 




4360 


4361 ± 9 ± 9 




(1s)t 3 - 


■T 4 


4418 


4421 ± 4 




(4s) 








4436 


(3d) 








4496 


(1«)t 4 - 


T 5 


4496 






(1«)t 5 - 


T 6 


4537 






(1*)t 6 - 


T 7 


4604 






(1«)t 7 - 


T S 


4673 


4664 ± 11 ±5 





TABLE IV: Calculated 1 charmonium masses, Meqm, 
from the energy dependent potential V(r) with a — 850 
MeV/fm, k = 100 MeV.fm and m c = 1348.5 MeV. Masses 
for experimental resonances, (M)pdg, have been taken from 
Masses from the Cornell potential, Mc or , are also shown 
for comparison. 
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jPC 

J 


Notation 


Meqm 


Mpdg 


M C or 






MeV 


MeV 


MeV 


0++ 


lp 


3456 


3414.75 ±0.31 


3456 


1++ 


lp 


3456 


3510.66 ±0.07 


3456 


2++ 


lp 


3456 


3556.20 ±0.09 


3456 


1++ 


2p 


3871.6 


3871.56 ±0.22 


3911 


0++ 


(lp)ri-Ob 


3898 












3915.5 ± 0.22 




2++ 


2p 


3903 


3929 ± 5 ± 2 


3911 


? : cr + 






3942!£ ± 6 




0++ 


(1p)To— To 








1++ 


(1p)t 1 -t 2 


4003 






1++ 


(1p)t 2 -t 3 








0++ 


(lp) T3 -T 4 


4140 












4143.0 ± 2.9 ± 1.2 










4156+Jg ± 15 




2 ++ 


(lp) Tl -T 2 


4140 






1++ 


(lp) T3 -T 4 


4190 






0,1,2++ 


(3p) 






4295 


0++ 


(lp) T4 -T 5 


4325 












4350.6tg;? ± 0.7 




TABLE V: 


Calculated 


J++ 


charmonium masses. 


Meq 



from the energy dependent potential V(r) with a = 850 
MeV/fm, k = 100 MeV.fm and m c = 1348.5 MeV. Masses for 
experimental resonances, (M) pdg, have been taken from [l3| . 
For p waves we quote separately the npo, npi and np2 states. 
Masses from the Cornell potential, Mc or , are also shown for 
comparison. The question mark indicates a plausible quan- 
tum number assignment for the experimental candidate (see 
Section III.B). 
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J 



PC 



Meqm Mgen 



[M Ti _ 1; M Ti ] MeV MeV 

^-i-Ti MeV (r 2 ) V2 (r 2 ) 1/2 () 



1"" 

To-Tx [0,4237] 1, 3046 3046 

(0.5) (0.5) 

3632 3632 

Is 1.2 

(0.9) (0.9) 

3743 3743 

Id 1.6 

(0.9) (0.9) 

4061 4066 

6s 12 

(1.2) (1.2) 

4137 4143 „ 

2d 15 

(1.3) (1.3) 

Ti - T 2 [4237, 43251 Is 4280 4290 32 

(1.5) (1.2) 

T 2 -T 3 [4325,43791 Is 4360 4374 43 

(1.8) 1.3 

T 3 -T 4 [4379,4470] Is 4418 4425 30 

1.5 1.3 

T 4 -T 5 [4470,4504] Is 4496 4513 51 

(2.3) (1.4) 

T 6 -T 6 [4504,4572] Is 4537 4545 35 

(1.7) (1.4) 

4604 4611 

T 6 -T 7 [4572,4648] Is 4DU4 4Dii 32 

(1.7) (1.5) 

T 7 -T 8 [4648,4685] Is 4673 4684 46 

(2.2) (1.5) 

TABLE VI: Calculated mixing angles 9(k) T . T . , generi- 

cally called 6, for the charmonium 1 __ EQM states. The 
corresponding interthreshold energy intervals, \MT i _ 1 ,MT i ~\, 
are shown. For the sake of completeness the EQM masses, 

Meqm, and the rms radii, (r 2 )^Q M , are also listed and 
compared to the generator masses, Mgen, and rms radii, 



GEN ' 
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PC Meqm Mgen 

[Mt^^MtA , MeV MeV 

Ti - 1_Ti MeV (r 2 ) 1/2 (r 2 ) 1/2 

(j = 1 ) \ / EQM V /GEN 



0++ 

To - Ti [0, 3730] lp 

Ti-T 2 [3730,3937] lp 34 

T 2 - T 3 [3937, 4014] lp 

T 3 -T 4 [4014,4224] lp ^™ " "' 26 

1.4 1.3 

T 4 -T s [4224,4405] lp 4325 4333 14 

1++ 

To - Ti [0, 3872] lp """" 3.5 



50 

(3.3) (1.1) 
Ti-T 2 [3872,4014] lp 40 ° 3 402 _ 9 46 

T 2 - T 3 [4014, 4080] lp 

T3-T4 [4080,4224] lp 37 

2++ 

T0-T1 [0,4014] lp """" """" 1.5 

(0-7) (0.7) 

3903 3911 
2p 17 
(1.2) (1.1) 

4140 4149 
Ti - T 2 [4014, 4224] lp ' 26 

(1.4) (1.3) 

TABLE VII: Calculated mixing angles &(k) T . T .i generi- 

cally called 6, for the charmonium J ++ EQM states. The 
corresponding interthreshold energy intervals, [Mt^,, MtJ, 
are shown. For the sake of completeness the EQM masses, 

Meqm, and the rms radii, (f 2 )^^- , are also listed and 
compared to the generator masses, Mgen, and rms radii, 



(fm) 


(fm) 


3456 


3456 


(0.7) 


(0.7) 


3898 


3916 


(1.3) 


(1.1) 




4083 




(1.3) 


4140 


4149 


1.4 


1.3 


4325 


4333 


1.5 


1.4 


3456 


3456 


(0.7) 


(0.7) 


3871.6 


3911 


(3.3) 


(1.1) 


4003 


4029 


(1.7) 


(1.2) 




4149 




(1.3) 


4190 


4206 


(1.6) 


(1.3) 


3456 


3456 


(0.7) 


(0.7) 


3903 


3911 


(1.2) 


(1.1) 


4140 


4149 


(1.4) 


(1.3) 



2\l/2 



GEN ' 



G, T PC\ ((cc)z=o + light-light meson) Mch 



Decay Channel MeV 



0~(1~~) J/ipfo(~ 500) ~ 3597 

J/V>/o(980) 4087 

V>(2s)/ (~ 500) ~4186 

J/^/ 2 (1270) 4372 

J/V>/i(1285) 4379 

J/V>/o(~ 1370) ~ 4467 

J/^/i(1420) 4523 

J/ -0/2 (1430) 4527 

</>(3s)/ (~ 500) 4540 

J/V>/o(1500) 4602 

J/V>/i(1510) 4615 

J/^/a(1525) 4627 

J/^/ 2 (1565) 4664 

V>(2s)/„(980) 4676 

0+(0,l,2++) J/iPuj(7S2) 3880 

J/^c/>(1020) 4117 

ij}(2a)u(782) 4469 



TABLE VIII: Two-body decay channels including I = hid- 
den charm states. For the calculation of their masses, Men, 
the quoted masses for (cc); = o and light-light mesons from [l3l ] 
have been used. When the mass of some light-light meson, 
/o(500) and /o(1370), is quite uncertain we have put a symbol 
of approximation (~) in front of its nominal value. 



Meqm Strong Decay Modes 
MeV 



r 



24 



Final Products 



3743 



(DD) P 
J/ipfo(~ 500) 



DD 

J/ipmv 



4061 



(D* D* , D S D S , 
DD*,DD)p 
J/ipfo(~ 500) 



(D*D*,D S D S 
DD* , DD) 
J/ipnir 



4137 



{DsD* a ,D*D*,D s Ds)p DsD* s ,D*D%D s Ds 
7/^/0(980) J/i>(irK,KK) 



{D* s D* S ,D S D* S ,D*D*) P D* D* , D S D* S , D* D* 

DD*tt 



(D£>i(2420), 
4280 DD!(2430))„ irt „ ai 
(V(2a)/o(~ 500), 
J/V/o(980)) 



DD*TT 

J/iP(tytt,KK)) 



D*D* (2400) D'Dtt 
4360 (V>(2s)/o(~ 500), (V>(2s)7T7r, 

(J/^/ 2 (1270))„ irtnai ) J/^tttt) 



4418 



L>*D5(2400) 
J/^(/o,i,a(~ 1300)) 



D* Dtt 

J/lp^AlT, 7T7r) 



4496 



L>*D*(2460) 
J/V(/o,i,2(~ 1300)) 



D* Dtt 
J /^{tt-k, 47r) 



D s A,i(2536) 
4537 (J/V(/o,i,2(~ 1400), 
</>(3s)/o(~ 500)) 



DsD'K 
(J/ip(KK~,nir), 



4604 



D s *D s i(2460) 
(J/^/o(1500), 
V(3a)/o(~ 500)) 



( J/ip(An, 7T7r) 
, 7/>(3s)7r7r) 



4673 



D:D s i(2536) 
V(2s)/ (980) 



D* 3 D*K 
V>(2s)(W,.Kif) 



TABLE IX: Dominant strong decays (subindex P for required 
P— wave) from the Quark- Antiquark and Threshold compo- 
nents for 1~~ EQM charmonium states with masses Meqm 
above the DD threshold. /o,i,a(~ 1300) and /o,i,a(~ 1400) 
stand for / (~ 1370) together with (/i(1285), / 2 (1270)) and 
(/i(1420), / 2 (1430)) respectively. The final products result 
from the decays of the second members of the decay channels 
(excepting D*'s and D^'s) as given in [l3| . 
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Meqm 
MeV 



Strong Decay Modes Final Products 



DD 

J/tpu}(782) 



DD 

J/i/;7r7r7r 



4140 



D*D* 
■7/^0(1020) 



D*D* 
J/iiKK 



4325 



J/V></>(1020) 
J/^w(782) 



D* a Dt 
J/ifjKK 



1++ 
3871.6 



(DD*) virtua i DDir 

(J/ljj0j(782)) v i r tual J/tpTTTTTV 



4003 



DD* 
J/V>w(782) 



DD* 

J/i/)7r7T7r 



4190 



D 3 D* 
J/^0(1O2O) 



D S D* 
J/ipKK 



2++ 
3903 



(DD)d 
J/^M 782) 



DD 



4140 



D*D* 
J/^0(1O2O) 



D* D* 
J/ipKK 



TABLE X: Dominant strong decays from the Quark- 
Antiquark and Threshold components for J ++ EQM charmo- 
nium states (subindex D for required D— wave) with masses 
Meqm above DD threshold. The final products result from 
the decays of the second members of the decay channels (ex- 
cepting D*'s and D*'s) as given in [l3| . 
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J PC Ti Mesonl - Meson2 (J^Jm-J Mt ( 

(0",1 + ) ?: 10978 



(i-,o+) 

?:H023 

(1",2+) 11072 
(0",1 + ) 11196 



0++ 



1++ 



2++ 



Ti 


B»B° 


(o- 


,0") 


10558 


T 2 


B*°W° 


(1- 




10650 


T 3 


B°B? 


(o- 




10734 


Ti 


B* a B^ 


(1- 


.1-) 


10830 


Ti 


B°W° 


(o- 


,1-) 


10604 


T 2 


B*°W° 


(1- 




10650 


T 3 


BgB* 


(o- 


.1-) 


10782 


T 4 


B* a B* 


(1- 


.1-) 


10830 


Ti 


B*°W° 


(1- 


.1-) 


10650 


T 2 


B* S B* 


(1- 


,1-) 


10830 



TABLE XI: Meson-meson thresholds for / = bb states. 



Ti 



S°Si(5721) 



■o 



B°B 1 {?) 

B*° B*(5732)° 
T 2 S*°B 1 (5721)° 

b*°bW)° 



T 3 B* U B 2 *(5747) 
Ta B°B 3 i(5830)° 



Notation 



Meqm 
MeV 



Mp DG 

MeV 
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Mcov 
MeV 



Is 


9459 


9460.30 ± 0.26 


9459 


2s 


10012 


10023.026 ± 0.31 


10012 


Id 


10157 


10163.7 ± 1.4 


10157 


3s 


10342 


10355.2 ± 0.5 


10342 


2d 


10438 




10438 


4s 


10608 


10579.4 ± 1.2 


10608 


3d 


10682 




10682 


5s 


10840 


10876 ± 11 


10841 


4d 


10898 




10902 


)t 1 -t 2 


10995 


11019 ±8 





(1s)t 2 -t 3 

(6s) 

(5d) 
(1s)t 3 -t 4 
(ld) T3 -T 4 
(2s)t 3 -t 4 



11039 



11090 
11130 
11140 



11053 
11105 



TABLE XII: Calculated 1 bottomonium masses, Meqm, 
from the energy dependent potential V(r) with a — 850 
MeV/fm, k = 100 MeV.fm and m b = 4793 MeV. Masses 
for experimental resonances, (M)pdg, have been taken from 
[ll|. Masses from the Cornell potential, Mc or , are also shown 
for comparison. 
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jPC 

J 


Notation 


Meqm 


Mpdg 


Mcor 








MeV 


MeV 


MeV 


0++ 


lp 




9920 


9859.44 ± 0.42 ± 0.31 


9920 


1++ 


lp 




9920 


9892.78 ± 0.26 ± 0.31 


9920 


2++ 


lp 




9920 


9912.21 ±0.26 ±0.31 


9920 


0++ 


2p 




10259 


10232.5 ±0.4 ±0.5 


10259 


1++ 


2p 




10259 


10255.46 ± 0.22 ± 0.50 


10259 


2++ 


■lp 




10259 


10268.65 ± 0.22 ± 0.50 


10259 


0++ 


3p 




10521 




10531 


1++ 


3p 




10526 


10530 ± 10 


10531 


2++ 


3p 




10528 




10531 


0++ 


0-P)t ± - 


-T 2 


10620 






0++ 


0)t 2 - 


-T 3 


10704 






1++ 


(1p)tx- 


-T 2 


_ 






1++ 


(1p)t 2 - 


-T3 


10708 






2++ 


(1p)tx- 


T 2 


10710 






0,1, 2+^ 


(4p) 








10768 


0++ 


(1p)ts- 


Ti 


10784 






2++ 


(2p)t x - 


T 2 


10815 






1++ 


(1p)t 3 - 


Ti 


10822 







TABLE XIII: Calculated J ++ bottomonium masses, Meqm, 
from the energy dependent potential V(r) with a — 850 
MeV/fm, k = 100 MeV.fm and m b = 4793 MeV. Masses 
for experimental resonances, (M)pdg, have been taken from 
[l3| . For p waves we quote separately the npo, npi and np2 
states. Masses from the Cornell potential, Mcor, are also 
shown for comparison. 
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Ti. 



Ti 



(i = l,...) 



[M Ti _ 1; M Ti ] 
MeV 





Meqm 


Mgen 




k 


MeV 


MeV 


on 


V / EQM 

(fin) 


V 1 GEN 

(fin) 


Is 


9459 
(0.2) 


9459 
(0.2) 





2s 


10012 


10012 





(0.5) 


(0.5) 


Id 


10157 


10157 





(0.6) 


(0.6) 


3s 


10342 


10342 





(0.8) 


(0.8) 


2d 


10438 


10438 





(0.8) 


(0.8) 


4s 


10608 
(1.0) 


10608 
(1.0) 


2.3 


3d 


10682 
(1.0) 


10682 
(1.0) 


2.2 


5s 


10840 
(1.2) 


10841 
(1.2) 


9.6 


4d 


10898 
(1.2) 


10902 
(1.2) 


14 


Is 


10995 
(1.2) 


10998 
(1.1) 


31 


Is 


11039 
(1.2) 


11042 
1.1 


29 


Is 


11090 
1.1 


11090 
1.1 


14 


Id 


11130 
1.5 


11133 
1.4 


23 


2s 


11140 
1.3 


11144 
1.2 


19 



T - Ti [0, 10978] 



T 1 -T 2 



T 2 -T 3 



T 3 -T 4 



[10978, 11023] 
[11023, 11072] 
[11072, 11196] 



TABLE XIV: Calculated mixing angles 6(k) T . l _ T . , generi- 

cally called 0, for the 1 EQM bottomonium states. The 
corresponding interthreshold energy intervals, [Mr i _ 1 , MtJ, 
are shown. For the sake of completeness the EQM masses 
M E qm and the rms radii (»' 2 )] 3 / 'q M are also listed and com- 
pared to the generator masses Mgen and rms radii (r 2 )^^ • 
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Meqm Mgen 



T;_ 



Ti 



MeV 



MeV 

/ EQM 

(fm) 


MeV 

V / GEP 

(fm) 


9920 


9920 


(0.4) 


(0.4) 


10259 


10259 


(0.7) 


(0.7) 


10521 


10531 


(1.0) 


(0.9) 


10620 


10628 


(1.1) 


(1.0) 


10704 


10711 


(1.1) 


(1.0) 


10784 


10789 


1.2 


1.1 


9920 


9920 


(0.4) 


(0.4) 


10259 


10259 


(0.7) 


(0.7) 


10526 


10531 


(1.0) 


(0.9) 




10669 




(1.0) 


10708 


10711 


(1.1) 


(1.0) 


10822 


10834 


(1.5) 


(1.1) 


9920 


9920 


(0.4) 


(0.4) 


10259 


10259 


(0.7) 


(0.7) 


10528 


10531 


(0.9) 


(0.9) 


10710 


10711 


(1.0) 


(1.0) 


10815 


10826 


(1.3) 


(1.0) 



0++ 
To - Ti 



Ti-T 2 
T2-T3 
T3-T4 

1++ 
To - Ti 



[0, 10558] lp 
2p 
3p 

[10558, 10650] lp 
[10650, 10734] lp 
[10734, 10830] lp 

[0, 10604] lp 
2p 



3p 

Ti - T 2 [10604, 10650] lp 
T 2 - T 3 [10650, 10782] lp 
Ti - T 4 [10782, 10830] lp 



2++ 

To - Ti [0, 10650] lp 
2p 
3p 

Ti - T 2 [10650, 10830] lp 
2p 



2.9 



24 



34 



34 



29 



2.5 



16 



22 



48 



11 



15 



32 



TABLE XV: Calculated mixing angles 0(k) T . T . , generically 
called 0, for the J ++ EQM bottomonium states. The corre- 
sponding interthreshold energy intervals, [Mr i _ 1 , Afr,], are 
shown. For the sake of completeness the EQM masses Meqm 
and the rms radii (^ 2 )£o M are also listed and compared to 



the generator masses Mgen and rms radii (r 2 ) 



2\l/2 



■G i 1 PC\ i(bb)l=o + light-light meson) M C h 



Decay Channel MeV 



_ (1 _ ~~) T(ls)/ (~ 500) ~ 9960 

T(ls)/„(980) 10450 

T(2s)/ (~ 500) ~ 10523 

T(ls)/ 2 (127Q) 10735 

T(ls)/i(1285) 10742 

T(ls)/o(~ 1370) ~ 10830 

T(3s)/ (~ 500) ~ 10855 

T(ls)/i(1420) 10886 

T(ls)/ 2 (1430) 10890 

T(ls)/ (1500) 10965 

T(ls)/i(1510) 10978 

T(ls)/a(1525) 10985 

T(2s)/ O (980) 11013 

T(ls)/ 2 (1565) 11022 

T(4s)/ (~ 500) 11079 

T(ls)/ a (1640) 11099 

0+(0, 1,2++) T(ls)w(782) 10243 

T(ls)0(lO2O) 10480 

T(2s)cj(782) 10806 

T(ls)w(~ 1420) 10880 

T(2s)<^(1020) 11043 

T(ls)^(1680) 11140 



TABLE XVI: Two-body decay channels including I = hid- 
den bottom states. For the calculation of their masses, Men , 
the quoted masses for (bb)i = o and light-light mesons from 
pH ] have been used. When the mass of some light-light me- 
son, /o(500), /o(1370) and w(1420) is quite uncertain we have 
put a symbol of approximation (~) in front of its nominal 
value. 



Mbqm Strong Decay Modes 
MeV 



Final Products 
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10608 



(BB) P 
(T(2a)/ (~ 500), 
T(ls)/ (980)) 



BB 

(T(2s)7T7r, 



[B*W, B~fF, BB) P 
10682 (T(2 S )/ (~ 500), 
T(ls)/ (980)) 



B*B~*, B~W, BB 

(T(2s)7T7r, 
r(la)(inr,KK)) 



10840 



(B a B*, B a B*,B a B a , 

B*W, BW) p 
(T(la)/o,i,a(~ 1300), 
T(3a)/o(~ 500)) 



(B a B a , B s B a , B s Bs 
B*W, B~B*) 
(T(ls)(7T7r,47r), 

T(3s)7T7r) 



10898 



{B* a B*,B a B*,B a B a , 

B*B*, BTF) p 
(T(la)/o,i,a(~ 1400), 
T(3a)/ (~ 500)) 



(BsB*, B a B a , B S B S , 

B*W, B~B*) 
(T(l«)(KKW,47r,7T7r), 

T(3s)7T7r) 



10995 



(BBl(5721, ?)) virtual 

(B*B*,B s B*,B a B s )p 
(T(la)/ 0> i,a(~ 1510), 
(T(2s)/ (980))„ lr t uai ) 



BB*n 
(B* B a , B a B* , B a B a ) 
{T{ls){KK, 4tt,tttv), 
T(2s)(tyty,KK)) 



11039 



(5*5^5732), 
BBi(5721,?)) 
(T(lfl)/ a (1565), 
T(2s)/ (980)) 



(B*B*7T, 
BB*7T) 
(T(ls)7T7r, 

T(2s)(7r7r,if3?)) 



11090 



B*B|(5747) 
(T(4a)/ (~ 500), 
T(ls)/ 2 (1640)) 



B*(B*ir, Btt) 

(T(4s)7T7T, 

T(ls)(ww,47r,ifi<')) 



11130 



B* 55(5747) 
T(ls)/ 2 (1640) 



B*(B*tt,Bit) 
T(1s)(uu,4it,KK) 



11140 



S* 55(5747) 
T(ls)/ 2 (1640) 



B*(B*tt,Btt) 
T{1s)(luuj,4tt,KK) 



TABLE XVII: Dominant strong decay channels (subindex 
P for required P— wave) from the Quark-Antiquark and 
Threshold components for 1 EQM bottomonium states 
with masses, Mbqm, above the BB threshold. Bi(5721, ?) 
stands forBi(5721) and Bi(?), /o,i,2(~ 1300,1400) for 
/o(~ 1370) with (/i (1285), f 2 (1270)') and (/i(1420), /a(1430)) 
respectively, and /o,i,a(~ 1510) for /o(~ 1500) with 
(/i (1510), / 2 (1525)). The final products result from the de- 
cays of the second members of the decay channels (excepting 
B*'s and B a 's) as given in [l3| . 
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Meqm 
MeV 



Strong Decay Modes Final Products 



0++ 
10620 

10704 

10784 

1++ 
10708 

10822 

2++ 
10710 

10815 



BB 
T(ls)0(lO2O) 
T(ls)w(782) 

B*B~* 
T(ls)^(1020) 

B s B a 
T(ls)<£(1020) 



B*B* 
T(ls)</>(1020) 

T(2s)o;(782) 
T(ls)0(lO2O) 



B*B* 
T(ls)0(lO2O) 

B*W 
T(2s)o;(782) 
T(ls)^(1020) 



BB 
T(ls)KK 
T(ls)7T7r7r 

B*B~* 
T(ls)KK 

B a B s 
T(ls)KK 



B*B* 
T(ls)KK 

B S B*, (B*B*) 
T(2s)7r7T7r 
T(ls)KK 



B*B* 
T(ls)KK 

B*TF 

T(2s)7r7T7r 
T{ls)KK 



TABLE XVIII: Dominant strong decay modes from the 
Quark- Antiquark and Threshold components for J ++ EQM 
bottomonium states with masses Meqm above BB thresh- 
old. The final products result from the decays of the second 
members of the decay channels (excepting B*'s and B^'s) as 
given in [l3| . 



